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We study light scalar mesons with particular focus on the ao(980) using 
lattice QCD with 2+1+1 dynamical quark flavors. To investigate the struc- 
ture of these scalar mesons and to identify, whether a sizeable tetraquark 
component is present, we use a large set of operators, including diquark- 
antidiquark, mesonic molecule and two-meson operators. We find that the 
low-lying states overlap essentially exclusively with two-meson states. This 
indicates that in the channels investigated no tightly bound four-quark 
states of either molecular or diquark-antidiquark type exist. 
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1. Introduction 

The nonet of light scalar mesons formed by o" = /o(500), k = ^"0(800), 
ao(980) and /o(980) is poorly understood. Compared to expectation all 
nine states are rather light and their ordering is inverted, which might in- 
dicate a strong tetraquark component. A detailed discussion of light scalar 
mesons can be found in [1]. They also have been discussed extensively on 
this conference (cf. e.g. the related publications [21 [3]). There are also var- 
ious tetraquark candidates among the heavy mesons, e.g. the rather light 
Z)*g(2317) and -Dsi(2460) mesons, whose masses seem to be difficult to re- 
produce theoretically using standard quark antiquark computations (cf. e.g. 

mum)- 

Here we report about the status of an ongoing long-term lattice QCD 
project with the aim to study possible tetraquark candidates from first 
principles. The focus of this talk is on the ao(980). Parts of this work have 
already been pubhshed [3 [U [9]. 

2. Lattice setup and four-quark creation operators 

We use gauge link configurations with 2-|-l-|-l dynamical quark flavors 
generated by the European Twisted Mass Collaboration (ETMC) [ini EH 
[T2| [T3l [T^ I15j . For the results shown in this talk disconnected diagrams 
have been ignored, which are technically rather challenging. An important 
physical consequence is that the quark number and the antiquark number 
are separately conserved for each flavor. Therefore, there is no mixing be- 
tween ■uu, dd and ss resulting in an r/^ meson with flavor structure ss instead 
of ?7 and 77' [8]. We are currently exploring efficient techniques to compute 
the relevant disconnected diagrams [9]. 

In the following we focus on the oo(980) sector, which has quantum 
numbers I{J^) = 1(0"^). As usual in lattice QCD we extract the low lying 
spectrum in that sector by studying the asymptotic exponential behavior of 
Euclidean correlation functions Cjk{t){{Oj{t))'^Ok{^)). Oj and Ok denote 
suitable creation operators, i.e. operators generating the oo(980) quantum 
numbers, when applied to the vacuum state. 

Assuming that the experimentally measured ao(980) with mass 
980 lb 20 MeV is a rather strongly bound four quark state, suitable creation 
operators to excite such a state are 

C^a'^fgso)" = E («"(^)75^(x)) (rf"(x)75Kx)) (1) 

X 

X 
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The first operator has the spin/color structure of a KK molecule (s(x)75ii(x) 
and (i(x)75s(x) correspond to a kaon K and an antikaon K at the same po- 
sition x). The second resembles a bound diquark antidiquark pair, where 
spin coupling via C75 corresponds to the lightest diquarks/antidiquarks (cf. 

e.g. m [HIE]). 

Further low lying states in this sector are the two particle states K + K 
and r]3 + TT. Suitable creation operators to resolve these states are 

<(9for'"'^'^ = (E^"(-)75-(x))(E^»75Ky)) (3) 

\ X y 

C^a^:Sor""'^''^ = (E-"(-)75Kx))(E^»75-(y)). (4) 

\ X y 

3. Numerical results an their interpretation 

We start by discussing numerical results for an ensemble with rather 
small spatial extension of L ~ 1.72 fm. This ensemble is particularly suited 
to distinguish two-particle states with relative momentum from states with 
two particles at rest and from possibly existing ao(980) tetraquark states 
(two-particle states with relative momentum have a rather large energy, 
because one quantum of momentum Pmm = 2t: /L 720 MeV). 

Figure [1^ shows effective mass plots from a 2 x 2 correlation matrix with 
a KK molecule operator ([1]) and a diquark-antidiquark operator ([2]). The 
corresponding two plateaus are around 1100 MeV and, therefore, consistent 
both with the expectation for possibly existing oo(980) tetraquark states and 
with two-particle K + K and r/s + vr states, where both particles are at rest 
{m{K + K)^ 2m{K) ^ 1198 MeV; m{rjs + vr) « m{r]s) + m{ir) « 1115 MeV 
in our lattice setup). 

Increasing this correlation matrix to 4 x 4 by adding two-particle K + K 
and 77s + vr operators (eqs. ^ and (jH)) yields the effective mass results 
shown in Figure [TJa. Two additional states are observed, whose plateaus are 
around 1500 MeV . . . 2000 MeV. From this 4x4 analysis we conclude the 
following: 

• We do not observe a third low-lying state around 1100 MeV, even 
though we provide operators, which are of tetraquark type as well as 
of two-particle type. This suggests that the two low-lying states are 
the expected two-particle K-\-K and ?7s + 7r states, while an additional 
stable ao(980) tetraquark state does not exist. 

• The effective masses of the two low-lying states are of much better 
quality in Figure [TJ) than in Figure [T^. We attribute this to the two- 
particle K + K and ry^ + vr operators, which presumably create larger 
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a) KK molecule , diquark (2x2 matrix) 



b) KK molecule , diquark , K+K two-particle , r\^+n {4x4 matrix) 
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Fig. 1. ao(980) sector, (L/a)'^ x (T/a) — 20^ x 48. a) Effective masses as functions 
of the temporal separation, 2x2 correlation matrix (operators: KK molecule, 
diquark-antidiquark, eqs. (3.2) and (3.3)). b) 4 x 4 correlation matrix (operators: 
KK molecule, diquark-antidiquark, two-particle K + K, two-particle rj^ + eqs. 
(3.2) to (3.5)). 



overlap to those states than the tetraquark operators. This in turn 
confirms the interpretation of the two observed low-lying states as 
two-particle states. 

• Analyzing the eigenvector components of the two low-lying states from 
Figure [TJ) we find that the lowest state is essentially exclusively of 
rjs + TT type, whereas the second lowest state is of -ftr -|- ^ type. On the 
other hand, the two tetraquark operators are irrelevant for resolving 
those states, i.e. they do not seem to contribute any important struc- 
ture, which is not already present in the two-particle operators. This 
gives additional strong support of the above interpretation of the two 
observed low lying states as two-particle states. 

• The energy of two-particle excitations with one relative quantum of 
momentum can be estimated according to 

m{l + 2,p = pmin) ~ ^m(l)2 +pI^^^ + ^m(2)2 + pl^^^ with 
Pmin = 27r/L. Inserting the meson masses corresponding to our lattice 
setup, m{K) Ki_ 599 MeV, m{r]s) 774 MeV and m(7r) 341 MeV, 
yields m{K + K,p = pmin) ~ 1873 MeV and 

m{r]s + TT,p = Pmin) ~ 1853 MeV. These numbers are consistent with 
the effective mass plateaus of the second and third excitation in Fig- 
ure [Ijj. Consequently, we also interpret them as two-particle states. 

We obtained qualitatively identical results, when varying the light quark 
mass and the spacetime volume [8j. 

Using exactly the same techniques, i.e. four-quark operators of tetraquark 
and of two-particle type, we also studied the K-sector (for details cf. [H]). 
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Again we find no sign of any four-quark bound state besides the expected 
two-particle spectrum (in this case K + tt states). Note that this result is 
in contradiction to a very similar recent lattice study of the k meson |18j . 
where an additional low lying four-quark bound state has been observed. 

4. Conclusions and future plans 

We have studied the ao(980) and the k channel by means of 2+1+1 fla- 
vor lattice QCD using four-quark operators of molecule, diquark and two- 
particle type. Besides the expected two-particle spectrum (two essentially 
non-interacting pseudoscalar mesons) no indication of any additional low 
lying state, in particular no sign of a four-quark bound state could be ob- 
served. This suggests that both the ao(980) and k meson have either no 
sizeable tetraquark component or they are rather weakly bound unstable 
states. To investigate the latter one needs to study the volume depen- 
dence of the two-particle spectrum in the corresponding sectors ( "Liischer's 
method", cf. e.g. [191 [201 [21]). Such computations are very challenging using 
lattice QCD, but first results have recently been published (cf. [22l[23]). We 
plan to perform similar computations with our setup in the near future. 
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